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The chirally symmetric overlap quark propagator is explored in Coulomb gauge for the first time. 
This gauge is especially well suited for studying the interrelation between confinement and chiral 
symmetry breaking, since confinement can be attributed to the infrared divergence of the vector 
dressing function of the inverse quark propagator. Using quenched gauge held configurations on 
a 20 4 lattice, the dressing functions of the quark propagator as well as the dynamical quark mass 
function are evaluated, also in the chiral limit. Chiral symmetry is artificially restored by removing 
the low mode contribution from the quark propagator. After removing enough low-lying modes, the 
dynamical quark mass function approaches the current quark mass in the whole momentum region 
and goes to zero in the chiral limit. However, the vector dressing function is unaffected by the 
low-mode removal. It follows that the quark energy dispersion is still infrared divergent and quarks 
in such a chiral symmetry restored phase are still confined within a hadron. 
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I. INTRODUCTION 

Spontaneous chiral symmetry breaking and confine¬ 
ment are the two main phenomena of low energy QCD. 
Chiral symmetry breaking goes in hand with dynamical 
mass generation, from which the low-lying hadron spec¬ 
trum can be understood. However, the underlying mech¬ 
anism of dynamical mass generation and the interplay 
between chiral symmetry breaking and confinement are 
still unresolved questions. 

The behavior of propagators allows to shed light on 
these unresolved issues. Coulomb gauge, where a confin¬ 
ing static quark potential arises naturally, is especially 
well suited for this purpose. From continuum mean-field 
studies using the so-called variational approach or Dyson- 
Schwinger equations it has been conjectured that confine¬ 
ment in this gauge reflects itself via infrared divergent 
quark propagator dressing functions, Refs. m, which, 
however, give an infrared constant dynamical mass func¬ 
tion, interpreted as the constituent mass of the quark. 
Only recently there has been a first attempt to trans¬ 
port this picture to full QCD using the lattice approach, 
Refs. [HE], giving hints that the quark self-energy is in¬ 
deed diverging in the low-momentum region. However, a 
quark propagator analysis with chirally symmetric lattice 
fermions is still lacking in Coulomb gauge. We attempt 
to close this gap. Overlap Dirac propagator studies in 
Landau gauge can be found in Refs. [6H9]. 

Continuum studies in Coulomb gauge are, in prin¬ 
ciple, appealing, for instance to explore the region of 
finite chemical potential m and to systematically ana¬ 
lyze the excited meson spectrum m However, although 
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some progress has been made in recent years m ep, 
continuum quark propagator studies in Coulomb gauge 
are still ansatz-dependent and far from being conclusive. 
Therefore, input from the lattice is highly needed. 

One of the main motivations of this work is to ex¬ 
plore, on the level of Green’s functions, what happens to 
the confinement properties of quarks, when chiral sym¬ 
metry is artificially restored by removing the low-lying 
eigenmodes from the Dirac operator in the valence quark 
sector. In recent years the effect of such an artificial 
symmetry restoration on the hadron spectrum has been 
analyzed, see Refs. [T41 - IT6] . From the exponential de¬ 
cay signal of hadron correlators it could indirectly be de¬ 
duced that the confinement properties stay intact. Here 
we want to give a more direct picture of the confine¬ 
ment issue by analyzing the low-mode truncated quark 
propagator in Coulomb gauge. The central question is 
how do the dressing functions change their shape when 
chiral symmetry is artifically restored and whether an IR- 
diverging dispersion relation is still possible. A first study 
of the low-mode truncated quark propagator has been 
made in Ref. m, however, where Landau gauge has been 
used, and a lattice fermion discretization, which only ap¬ 
proximately fulfills the Ginsparg-Wilson equation. 

Having a chirally symmetric lattice Dirac operator at 
hand, is central for our purpose. We use the overlap 
Dirac operator, Refs. mum, which is an explicit solu¬ 
tion to the Ginsparg-Wilson equation. Besides its nice 
properties, like the occurence of exact zero modes, it ad¬ 
mits a clear prescription to extract the dressing functions 
of the quark propagator. After identifying the overlap 
lattice momenta, no further tree-level corrections or im¬ 
provements have to be performed, which could lead to 
ambiguities in the result. On the other hand, since over¬ 
lap fermions are numerically very costly, we are restricted 
to rather small lattices. Nevertheless, we can enter the 
small momentum region of the dynamical mass function 
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and evaluate the low energy chiral properties of QCD. 

The organization of the article is as follows: in 
Chapter [Tl] the overlap Dirac operator is defined, the free 
case is discussed and the overlap lattice momenta are de¬ 
rived. In Chapter |III| the lattice setup, the numerical 
implementation of the overlap Dirac operator and the 
gauge fixing procedure are described. In Chapter IV the 
Coulomb gauge quark propagator is introduced and the 
numerical results for the dressing functions are presented. 
In Chapter [V] the Dirac low mode removal is discussed, 
its relation to confinement is clarified and the dressing 
functions for different truncation steps are shown. In 
Chapter |VI| the main results of this work are summar¬ 
ized and a short conclusion and outlook are given. 



II. OVERLAP DIRAC OPERATOR AND FREE 
PROPAGATOR 


Figure 1. Overlap lattice momentum q(p) for six different 
masses in units of MeV (a = 0.2 fm). 


The massless overlap Dirac operator is defined as 

D( 0) = p (1 + 75 sign [H w (-p)]) , (1) 

with Hw(—p) = a^D^^—p) the (dimensionless) Her- 
mitian Wilson-Dirac operator and p a negative mass 
parameter. The operator 0 fullfills the Ginsparg-Wilson 
equation [ 2 U] 

D( 0)75 + 75-0(0) = -D(0) 7 5-D(0) , (2) 

P 

and therefore represents exact chiral symmetry on the 
lattice, see Refs. [2TH25] . The eigenvalues of D(0) lie on 
a circle in the complex plane with radius p. The free 
massless overlap quark propagator in momentum space, 
denoted as S^{p), can be decomposed as 

S {0) {p) = -i'Yp.C^p) + Y p > ( 3 ) 


with 



and the momenta k^ = sin(p^a),^ = 2sin(p M a/2). In 
contrast to its continuum counterpart, the free massless 
quark propagator has an additional term 1 /( 2 p), which 
enters with the unit matrix in Dirac space. It is allowed 
in the Ginsparg-Wilson equation ([2|. 

According to Ref. 0, we redefine the massless overlap 
quark propagator as 


S = S- 


1 

2 P ’ 


( 5 ) 


which corresponds to imposing continuum chiral sym¬ 
metry for overlap fermions: {S, 75 } = 0. The free mass¬ 
less lattice quark propagator then has the same Dirac 


structure as the free continuum propagator. No addi¬ 
tional improvement techniques have to be imposed, which 
could introduce a small arbitrariness, depending on the 
particular procedure, see Ref. [BJ. This is a great advant¬ 
age of choosing overlap fermions. 

The massive overlap Dirac operator follows from 
Eq. (jTJ) via 

D(m 0 ) = ^1 — D(0) + mo , (6) 

with to 0 the quark mass parameter. From the massive 
free (inverse) overlap quark propagator, which has the 
structure 

(S (0) ) (p) = + lm > ( 7 ) 

we identify the overlap lattice momenta q M and the cur¬ 
rent quark mass m 

4 p 2 ( \J k l + A 2 -f- A) _ too 

g »-(2p-m 0 ) kl 1 

( 8 ) 

As a first test of our code we compare the numerically 
evaluated free overlap quark propagator with the ana¬ 
lytic formulas (|8]) , see Figs, [lj [2] A similar plot of the 
overlap lattice momenta q^ip) can be found in Ref. [B]. 
Throughout this work we perform a cylinder cut (Z¥j on 
all data, that is we evaluate the numerical functions only 
within one unit of spatial momentum from the diagonal 
of four-momentum space. 


III. SETUP OF THE CALCULATION 

A. Gauge field configurations 

We use quenched Liischer-Weisz [25] gauge field config¬ 
urations which we generated with the Chroma software 
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in which the linear definition of the gauge fields is used, 


A i(x) 


Ui(x ) - Ut(x )+ 


2iag 0 


traceless 


( 12 ) 


where g q ist the coupling constant. 

We maximize F g [U] with the cuLGT code [31] and the 
overrelaxation algorithm [32] . which in the case of Cou¬ 
lomb gauge can be applied to each time-slice independ¬ 
ently. 

A measure of the quality of the gauge fixing is the 
average Z^-norm of the gauge fixing violation A 9 ^ 0 

0(z 4 ) = tr [ A9 ( a: > a; 4)A 9 (a;, x 4 ) f ] , (13) 

-' v S Jv c x 


Figure 2. Current quark mass m for free overlap fermions. 


package [26] and QDP-JIT [27, [25] • We adopted an in¬ 
verse coupling value /3± = 7.552 which corresponds to a 
lattice spacing a = 0.2 frn [25] The size of the lattice is 
Ns x Nt = 20 1 * 3 * * x 20 and the value of the average plaquette 
is 0.5767. 

We use an ensemble of 96 well separated configurations 
and evaluate the quark propagator for six current quark 
masses chosen at m = 85,92,99,113,137,173 MeV. The 
quark propagators are obtained by inverting the over¬ 
lap Dirac matrix for each configuration on point-sources. 
Subsequently, the quark propagators are transformed to 
momentum space after which the extraction of the dress¬ 
ing functions is performed according to the description 

in mesa. 


where the sum runs over all spatial sites x within one 
time-slice and Ns is the number of spatial lattice sites. 
We have chosen 0 < 10~ 12 as the stopping criterion for 
the overrelaxation algorithm. 


In Coulomb gauge, maximizing Eq. (10) leaves the 
temporal links U±{x) unfixed, i.e., a residual gauge free¬ 
dom with respect to space independent gauge transform¬ 
ations g(x±) is left. One possible choice to fix the residual 
gauge in the continuum is to require 


c>4 J d 3 x A 4 (x) = 0 . 


(14) 


We will use the lattice version of Eq. (14), the so-called 
integrated Polyakov gauge if 551 . 


C. Overlap Dirac operator 


B. Gauge fixing 


The continuum Coulomb gauge condition, 

d l A i = 0, * = 1,2,3, (9) 

can be realized on the lattice by maximizing the gauge 
functional 

F g [U] =5R^tr [g(x)Ui(x) g(x + tf] (10) 

i,x 


with respect to gauge transformations g(x) £ SU(3). The 
gauge links maximizing Eq. (10) satisfy the discretized 
Coulomb gauge condition 


A9 (a) = ( A iW ~ <>-*)) = o, ( n ) 

i 


1 Both, the Liischer-Weisz gauge action and the overlap fermi¬ 

onic action are improved, to orders 0(a?) and 0(a), respectively, 

which allow the usage of a rather large lattice spacing. Moreover, 

this study focuses on the infrared behaviour of QCD, where the 

ultraviolet cutoff is of minor importance. 


The (valence) overlap Dirac operator is implemented 
in the following way: the low-lying part of sign[ZJ\v] 
is calculated via the spectral decomposition. For the 
high-lying part the Chebyshev approximation is used, see 
Ref. [33]. Eigenvalues and -vectors are extracted with the 
ARPACK routines, Ref. [[35]. For the inversion of the 
overlap Dirac operator the biconjugate gradient stabil¬ 
ized multimass solver is used, see Ref. [36]. The Wilson- 
Dirac mass parameter p is set to 1.6. 

In order to reduce the cost for the overlap Dirac op¬ 
erator, three sweeps of stout smearing m are applied 
on the gauge field configurations. This procedure is es¬ 
pecially important when computing the low-lying eigen- 
modes of D. After three sweeps of stout smearing the 
configurations are smooth enough so that the low modes 
of the Hermitian Wilson-Dirac operator are suppressed, 
see Ref. [35] for details. This results in a speed-up of the 
overlap construction by a factor of three. 

The propagators are computed on unsmeared and 
stout smeared configurations. The low-mode truncated 
propagators, i.e. where the low-lying eigenvalues are re¬ 
moved from the propagators, are computed on smeared 
configurations only, since for the evaluation of the lowest 
overlap eigenvalues smearing is essential to reduce the 
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computational cost. In Chapter |IV C we analyze the ef¬ 
fect of the smearing procedure on the quark propagator 
dressing functions. 


IV. NON-PERTURBATIVE QUARK 
PROPAGATOR IN COULOMB GAUGE 

A. Introduction 


In Coulomb gauge, due to the breaking of Lorentz in¬ 
variance, four independent scalar dressing functions for 
the (Euclidean) quark propagator occur: 

S~ 1 (p 1 p i ) = i'YiPiAsip ,p 4 ) + *74P4 j4t(p, Pa) (15) 

+'HP4'YiPiA D (p,p i ) + lB(p,p4) , 


with A s ,At,A 0 ,B referring to spatial , temporal , mixed 
and scalar dressing functions i39:, respectively. Here p 
denotes 3-momentum. For free fermions in the (Euc¬ 
lidean) continuum the dressing functions are = A T = 
1 ,B = mo and A B = 0. 

From a recent study in (lattice) Coulomb gauge, 
Ref. [4], it has been shown that the dressing functions 
in Eq. ( fl5| ) are independent of the fq-component, which 
is also supported by our results with overlap fermi¬ 
ons. Moreover, we also find that the mixed component 
A 0 (p,p 4 ) vanishes. For the smallest momenta, however, 
the data suffer from large statistical fluctuations. To in¬ 
crease statistics, we average the dressing functions over 
the lattice cubic symmetries, the parity symmetry of the 
quark propagator and finally over the jq-component. 

From mean-field studies in continuum Coulomb gauge 
it is suspected that, due to the presence of the non- 
Abelian color-Coulomb potential, the scalar and vector 
dressing functions diverge for \p\ —> 0, see Refs. @131- 
The dynamical quark mass M(p), which is defined as 


M{p) 


B (P) 

Mp) ’ 


(16) 


however, reaches a constant value for \p\ —> 0, identified 
as the constituent quark mass. The infrared divergencies 
of B(p) and A s (p) cancel each other, so that M(p) is 
infrared finite, see Refs. mm for details. It is this di¬ 
vergence property of scalar and vector dressing functions 
that makes Coulomb gauge special, since it can be re¬ 
lated to the confinement properties of quarks, as we will 
show in the next section. 

We stress that in the continuum non-linear integral 
equations (Dyson-Schwinger equations) an infrared reg¬ 
ulator has to be introduced to make the integrals, due to 
the non-Abelian color-Coulomb potential, well-defined. 
The divergence of the dressing functions appears for van¬ 
ishing infrared regulator. For a finite lattice the dress¬ 
ing functions then should approach a constant value for 

IpI = °- 

On the lattice, the regularized quark propagator is ex¬ 
tracted which depends on the lattice cutoff a. The reg¬ 
ularized quark propagator Sl [p] a) can be renormalized 


at the renormalization point £ with the momentum in¬ 
dependent quark wave-function renormalization constant 
Z 2 U;a), 

S L {p-a) = Z 2 {^a)S^ 1 p). (17) 

While the mass function M(p ) is independent of the 
renormalization point £, the wave-function renormaliz¬ 
ation function 


Z*(p) = 1 Ms(p) ( 18 ) 

differs at different scales but can be related by multiplic¬ 
ation with a constant. 


B. Numerical results 


Here we discuss our results for the dressing functions 
of the inverse quark propagator S 1-1 ^), Eq. (15) as well 
as the dynamical quark mass M(p ), Eq. (16), and the 
wave-function renormalization function Z^(p), Eq. (18). 

Most interesting are the spatial and scalar dressing 
functions, A s (p) and B(p). For large momenta A s (p) 
goes to unity and B(p) approaches the current quark 
mass consistent with asymptotic freedom. Around one 
GeV both dressing functions increase for small momenta, 
which goes in hand with spontanteous symmetry break¬ 
ing and dynamical mass generation, see Figs. [3j [4j We 
also show a simple linear extrapolation to the chiral limit. 
For massless quarks B(p) goes to zero for \p\ —► oo. We 
note that for zero momentum B(p) reaches a finite value 
which is most likely a finite size effect. The vector dress¬ 
ing function A s (p) is not accessible at zero momentum. 
For the scalar dressing function B(p) larger current quark 
masses lead to larger values for small momenta. The 
spatial dressing function A s (p) only has a mild mass de¬ 
pendence. All values for the different masses lie within 
error bars. In Fig. [5] the renormalization function Z^{p), 
renormalized at £ = 6 GeV, is shown. To verify that 
A s (p) is indeed IR diverging, or alternatively that Z^(p') 
goes to zero for \p\ —»• 0, the infinite volume limit has to 
be taken. This is left to a further investigation. 

The dynamical quark mass function M(p) is shown 
in Fig. [6] For large momenta the mass function should 
approach the current quark mass m, which is, due to 
the rather coarse lattice, only approximately seen in our 
results. In the chiral limit it goes to zero, which is, 
within error bars, also approximately seen in our res¬ 
ults. Between one and two GeV the function starts to 
increase and acquires large values for small momenta. 
The smaller the current quark mass, the more dynamical 
quark mass is generated. This means that for larger cur¬ 
rent quark masses the amount of IR mass generation is 
smaller. This becomes obvious by comparing largest and 
smallest current quark masses m = 173 MeV and m = 85 
MeV, respectively. At the smallest accessible momentum 
their values are split only by around 50 MeV. This beha¬ 
viour reflects itself in the strong infrared dynamical mass 
generation for chiral quarks. 
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Figure 3. Spatial component A s (p) for several quark masses Figure 5. Z^(p) renormalized at £ = 6 GeV for smallest quark 

and in the chiral limit. mass m = 0.085 GeV. 



Figure 4. Scalar component B(p) for several quark masses 
and in the chiral limit. 



Figure 6. Dynamical mass function M(p ) for several quark 
masses and in the chiral limit. 


Due to the relatively large lattice we reach IR momenta 
down to around 300 MeV which is, compared to former 
Coulomb gauge studies, quite an improvement. At \q\ « 
300 MeV the dynamical quark mass is already at 210 
MeV for massless quarks. It is thus very likely that a 
constituent quark mass around 300 MeV is reached. 

We also note that the dynamical quark mass M{p) 
should only be affected by vacuum fermion loops (dy¬ 
namical configurations) at the percentage level, Ref. [3]. 
Hence, when using costly dynamical overlap configura¬ 
tions we do not expect much difference for the dynamical 
mass function. 

The mixed dressing function A D ( p ) seems to vanish for 
all current quark masses, see Fig. [7] For small momenta 
the error bars are too large to make a precise statement. 

Finally, we comment on the temporal dressing func¬ 
tion A T (p), which is only well-defined when the residual 
gauge invariance with respect to space independent gauge 
transformations is fixed. However, it is not influenced by 
the Coulomb gauge condition. The additional gauge free¬ 


dom is fixed via the Integrated Polyakov gauge. In this 
gauge A T (p) goes to small non-zero values for large mo¬ 
menta. For small momenta the error bars are too large 
to make a precise statement, see Fig. [8j However, a di¬ 
vergent behavior, as for the dressing functions fixed by 
Coulomb gauge, does not seem to take place. We note 
that the qualitative behavior is consistent with the find¬ 
ings in Ref. [3]. For large and intermediate momenta no 
mass dependence is seen. For small momenta there is 
a tendency that smaller masses approach larger values. 
However, one has to be careful, since the error bars are 
large in this region to make a precise statement. The 
temporal dressing function should also be clarified in the 
continuum approach, which, until now, has not been con¬ 
sidered. 

We stress that for the confinement properties the dress¬ 
ing function A s (p) is the quantity of central interest. This 
issue will be clarified in the next section. 
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Figure 7. Mixed dressing function A D (p) for several current 
quark masses. 



Figure 8. Temporal dressing function A T (p) for several cur¬ 
rent quark masses and in the chiral limit. 


C. Comparison with propagator on stout-smeared 
configurations 

In order to understand the role of smearing for the 
quark propagator in Coulomb gauge, we compare the 
dressing functions on stout-smeared configurations with 
the unsmeared case. We use the same mass paramet¬ 
ers m o in both cases. We do not interpolate between 
the mass values, as done in Ref. [J2], where the effect 
of smearing on the Landau gauge quark propagator was 
under consideration. 

All dressing functions show the same qualitative beha¬ 
vior as in the unsmeared case. The functions A s (p), B(p) 
and A T (p) still increase for small momenta, see Figs. |9j 
mm However, the dressing functions approach dif¬ 
ferent low-momentum values as in the unsmeared case. 
The functions A s (p) and B{p) reach smaller values for 
low-momenta, whereas A T (p) reaches larger values. In 
the mid-momentum regime A s (p) gives smaller values 
as compared to the unsmeared case, whereas the func- 


Figure 9. Comparison of spatial component A s (p) with three 
levels of stout smeared and unsmeared configurations for the 
mass value m = 99 MeV. 


tion B{p) gives nearly the same values as in the un¬ 
smeared case. The scalar dressing functions start to devi¬ 
ate from each other around 1 GeV. For A T (p) the values 
are slightly larger in the whole momentum regime. 


We also observe, that the error bars for the IR values 
for all dressing functions are drastically reduced com¬ 
pared to the unsmeared case. This is especially import¬ 
ant for the temporal dressing function A T (jp), see Fig. 11 
which gives now a clearer IR behavior, as for the un¬ 
smeared case, see Fig. [8] for comparison. From this result 
we can conclude that the mass dependence of the tem¬ 
poral dressing function A T (p) is small. The fact, that the 
error bars are reduced for the IR values, seems to be a 
direct consequence of the smoothing procedure. 


In Fig. 12 we compare M(p) for smeared and un¬ 
smeared configurations for in = 99 MeV and in the chiral 
limit. In contrast to the dressing functions, the dynam¬ 
ical quark mass M(p) agrees quantitatively, within the 
error bars, for smeared and unsmeared configurations. In 
the smeared case it reaches slightly larger infrared values 
for finite current quark mass m, which is due to the fact, 
that A s (p) in the IR is smaller than for unsmeared config¬ 
urations. Clearly, as pointed out in Ref. [42], the overlap 
mass parameter mo on the smeared configurations gives a 
different renormalized quark mass. However, in the chiral 
limit smeared and unsmeared propagators are nearly in¬ 
distinguishable, and a constituent mass of around 300 
MeV is reached. This is a remarkable result: although 
scalar and vector dressing functions give quite different 
values for the smeared propagator, the dynamical quark 
mass M(p) is not really affected by the smearing pro¬ 
cedure. The differences cancel each other in the ratio 


Eq. (16). 
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Figure 10. Comparison of scalar component B(p) with three 
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Figure 11. Temporal dressing function A T (p) with three levels 
of stout smeared and unsmeared configurations for the mass 
value m = 99 MeV. 
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V. CONFINEMENT IN COULOMB GAUGE 
AND ARTIFICIAL CHIRAL SYMMETRY 
RESTORATION 


Figure 12. Dynamical mass function M(p) with three levels 
of stout smeared and unsmeared configurations for the mass 
value m = 99 MeV (open symbols) and in the chiral limit 
(full symbols). 


integration for the interacting quark propagator S(p,p 4 ), 
Eq. (15), yields 


S(p) = 


B(p) — ■ pA s (p) 

Mp) 


( 20 ) 


with the dispersion relation of a confined quark 

U}(p) = A T (p)yJp 2 A 2 (p) + B 2 (p) . (21) 

If scalar and vector dressing functions, B(p) and A$(p ), 
respectively, are IR divergent, the energy dispersion to(p) 
is IR divergent as well. This gives rise to the physical 
picture of quark confinement: the excitation energy of a 
single quark is infinite. It has to be stressed that such 
a mechanism of quark confinement via an IR divergent 
quark dispersion relation is special for Coulomb gauge 
and absent, for instance, in Landau gauge. 


Using the dynamical quark mass M(p), Eq. (16), we 
can rewrite Eq. © as 

w(p) = A T (p)A s (p)^/p 2 + M 2 (p) . (22) 


We clarify how confinement is related to the infrared 
divergent quark dressing functions in Coulomb gauge. 
Such a connection has first been discovered in a con¬ 
tinuum Coulomb gauge model, Refs. m- Only recently 
it has been explored in lattice calculations as well, see 
Refs. [H[5]. 

Integrating the free (Euclidean) quark propagator over 
the component p 4 gives the free static quark propagator 

dp 4 1 too — i-f ■ p 

27r ij ■ p + iy 4 P 4 + wo 2w(p) 

(19) 

where in the denominator the free-particle dispersion re¬ 
lation u>(p) = p 2 + to 2 appears. Performing the P 4 - 


Supposing that M(p) and At(p) are constant for |p| —>■ 
0, the IR-divergence of u>(p) is fully included in the in¬ 
frared divergence of As(p). Hence, the vector dressing 
function Ag(p) is crucial to link the confinement proper¬ 
ties of quarks to chiral symmetry breaking. In Ref. [3] 
first indications of an IR-diverging quark dispersion rela¬ 
tion uj(p) on the lattice have been observed. 

Now an interesting question arises: if removing the 
chiral condensate from the quark propagator, is the en¬ 
ergy dispersion w(p), Eq. ( plj ), still infrared divergent? If 
yes, then confinement is intact, although chiral symmetry 
is artificially restored in the vacuum. We will analyze this 
question numerically. Via the Banks-Casher relation [43] 
the chiral condensate is connected to the low-lying eigen¬ 
values of the Dirac operator. These low-lying eigenvalues 







































are removed via the prescription 


ok 

^TRUNC 


= S V 


Y 


(23) 


with Strung an d S FULL denoting low-mode truncated and 
full propagators, respectively. Here A* are the eigenval¬ 
ues, | Vi) the eigenvectors of the overlap Dirac operator 
and k denotes the number of truncated modes. For in¬ 
stance, k = 16 means subtracting the 16 lowest modes 
from the quark propagator. For our 20 4 lattice we use 
truncation steps from 16 to 256, corresponding to remov¬ 
ing all eigenvalues of the massless Dirac operator spec¬ 
trum up to values of 6 — 138 MeV, see Tab. [T] for the 
single steps. After removing enough modes, the quark 
condensate is zero. The truncated propagators Sprung 
are evaluated for all six mass values and the chiral limit 
is considered as well. 


k 

cutoff [MeV] 

16 

6 ± 3 

32 

14 ±3 

64 

29 ±3 

96 

46 ±4 

128 

63 ±5 

256 

138 ±6 


Table I. The truncation steps k and the corresponding ei¬ 
genvalue cutoffs of the massless Dirac operator at which we 
evaluate the quark propagator. The cutoff values are averages 
over all configurations and the errors are standard deviations. 


It is important to note that while the truncation 
leaves unitarity intact, by construction it violates the loc¬ 
ality of the Dirac operator to some extent. In [H] it has 
been shown, however, that the violation of locality of 
a truncated Wilson type Dirac operator is many orders 
of magnitude smaller than the ultra local contribution. 
Moreover, in Ref. US] it has been demonstrated that, at 
finite lattice spacing, the (untruncated) overlap operator 
reveals nonlocal contributions which fall of exponentially 
as a function of the distance. Therefore we are confident 
that the violation of locality resulting from the low mode 
truncation Eq. (23) does not affect our results. 

After chiral symmetry restoration, the quark condens¬ 
ate (ipip) and therefore the dynamical quark mass M (p) 
and scalar dressing function B(p) vanish in the chiral 
limit. The quark dispersion relation oj(p ), Eq. (22), then 
reads 


Wtrunc (p) = A T (p)A s (p) |p| . (24) 

The important question now is, what happens with the 
vector dressing function Ag (p) after chiral symmetry res¬ 
toration. This goes in hand with the question, whether 
confinement survives artificial chiral symmetry restora¬ 
tion 


2 It can be speculated that Ag(p) serves an an order parameter 


We note that in Ref. [H], within a continuum Coulomb 
gauge model, it is demonstrated that with removing all 
the momenta relevant for the chiral condensate, confine¬ 
ment is still intact. It is shown that the vector dressing 
function Ag(p) is still infrared divergent, while the scalar 
dressing function B(p) goes to zero. In the next section 
we will explore numerically, if it is possible to generalize 
this picture to full QCD. 

In the remainder of this section we clarify how the 
dynamical quark mass M (p), a gauge variant quantity, 
is connected to the chiral condensate (ip(x)ip(x )) which, 
in contrast, is gauge invariant. 

Using Eq. (20) the chiral condensate, which is given as 


the trace of the static quark propagator, can be expressed 
as 


(ip ( x)ip a (x )) = —2 Ac 


M (p) 


d 3 p 

(27t) 3 A t (p)Yp 2 * * 


M 2 (p) 
(25) 


We note that the integral on the right-hand side is only 
well-defined in the chiral limit. From this formula it fol¬ 
lows that the chiral condensate vanishes when M(p) = 0, 
assuming that A T (p) stays finite. 

The expression ( [25] ) could be used to relate Dirac ei¬ 
genvalues, used on the lattice to artifically remove the 
chiral condensate, and quark momenta, used in func¬ 
tional methods like Dyson-Schwinger equations for this 
purpose. For instance, it is still not fully understood why, 
after only a small amount of low-modes are removed, 
where the chiral condensate is not yet zero, mesons 
already appear in chiral multiplets. New insights on 
this issue could be gained by comparing these two ap¬ 
proached 


A. Numerical results 


For each current quark mass to we evaluate the dress¬ 
ing functions for all six truncation steps. The dressing 
functions are shown in the chiral limit. The mass func¬ 
tion M(p) and scalar dressing function B(p) are also 
shown for the current quark mass to = 99 MeV. Our 
conclusions apply for all mass values. 

We start with discussing the dynamical quark mass 
M(p), Fig. 13 As expected from chiral symmetry consid¬ 
erations, it decreases for increasing number of truncated 


of confinement in Coulomb gauge. This issue will be explored 

further in future studies, for instance, by studying the quark 
propagator in the deconfinement region. 

3 The procedure would be as follows: after removing a small 
amount of low-lying modes, the constituent quark mass M (0) has 
a certain non-zero value below the QCD value around 300 MeV. 
In a functional approach the integration momenta |p| < |p| cr.it in 
Eq ,< |25| ) should be excluded and |p|crit adjusted to get the same 
value for M( 0). Then the meson spectrum should be evaluated 
by solving the Bethe-Salpeter equation. 
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Figure 13. Dynamical mass function M ( p ) for the untrun¬ 
cated case and for six truncation steps for the current quark 
mass m = 99 MeV. 



Iql [GeV] 

Figure 14. Dynamical mass function M(p) for the untrun¬ 
cated case and for six truncation steps in the chiral limit. 


Dirac eigenmodes. Already after the first truncation step 
(k = 16), M(p) has lost some IR strength. After remov¬ 
ing k = 256 modes the function appears close to the 
current quark mass, M(p) ss to. In the chiral limit this 
results in M(p) approaching a small (but non-zero) IR 
value, see Fig. [l4j after k = 256 low modes have been re¬ 
moved. Correspondingly, the chiral condensate is small 
but not yet zero at that stage. 


The spatial dressing function As(p), however, does not 


decrease after removing the low-lying modes, see Fig. 15 


Instead, the small momentum values even increase after 
removing the low-lying eigenvalues. While for large and 
intermediate momenta (around 1 GeV) the values for 
each truncation step lie on top of each other. The former 
is one of the crucial observations of this study. It is a clear 
indication that an infrared divergent dispersion relation 
w(p), Eq. (24), survives the artificial chiral symmetry 


Figure 15. Spatial dressing function As(p) for the untrun¬ 
cated case and for six truncation steps in the chiral limit. 
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Figure 16. Scalar dressing function B(p) for the untruncated 
case and for six truncation steps for the current quark mass 
m = 99 MeV. 


restoration 0 

The scalar dressing function B(p) decreases under low¬ 
mode truncation, see Fig. |16| consistent with chiral sym¬ 
metry considerations. However, for non-vanishing quark 
mass B(p) increases also at the largest truncation level 
k = 256 quite strongly towards the IR. This behavior 
can be understood by considering the mean-field Dyson- 
Schwinger equations of the truncated theory (Ref. [46] 1: 
the bare quarks still interact with gluons via the color- 
Coulomb potential. Only in the chiral limit B{p) van¬ 
ishes after removing enough low inodes, since the Dirac 


4 In a next step the IR behavior of the spatial dressing function 
should be analyzed in detail, for instance, it should be checked, 
if it obeys a power-law Ag(p —» 0) ~ |p| — “, and understood, 
why after truncation, the value for a becomes even larger. For 
such an analysis the continuum Coulomb gauge model could be 
well-suited. 
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Figure 17. Scalar dressing function B(jp ) for the untruncated 
case and for six truncation steps in the chiral limit. 
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Figure 18. Temporal dressing function At(p ) for the untrun¬ 
cated case and for six truncation steps in the chiral limit. 


structure 1 is then not allowed. In our case, after remov¬ 
ing k = 256 modes, B(p) has lost more than the half of 
its IR strength, see Fig. [IT] But it is still non-zero, which 
again confirms that the chiral condensate is non-zero also 
at the largest truncation levelQ 

The temporal dressing function At(p ) is also allowed 
in the chirally symmetric scenario. It stays non-zero after 
truncation as shown in Fig. |18| It behaves similar to the 
spatial dressing function fig (p). It does not change under 
truncation for low and intermediate momenta. For small 
momenta it seems to increase, although the gauge noise 
is too large to make a conclusive statement. 


5 In Coulomb gauge the infrared divergencies, which occur for spa¬ 
tial and scalar dressing functions, cancel each other in the dy¬ 
namical mass function and the observable quantites. Since in the 
chirally symmetric phase the dynamical quark mass M (p) wants 
to approach the current quark mass m in the IR, B (p) and A s (p) 
have to change their IR behavior, so that M (p) is still finite. 


The mixed component A d ( p ), which does not occur in 
the untruncated case, is still vanishing for the truncated 
case and therefore not shown here. 

Let us underline the main outcome of this numerical 
study. After removing the low-lying eigenmodes from 
the quark propator, which are responsible for the form¬ 
ation of the chiral condensate, the Lorentz-vector dress¬ 
ing function, which enters the quark dispersion relation 
w(p), Eq. (24), stays intact and even becomes stronger 
divergent in the IR. This gives the Coulomb gauge ex¬ 
planation of confinement after artifical chiral symmetry 
restoration: due to the persistence of the vector dress¬ 
ing function Ag(p), the quark dispersion relation w(p), 
Eq. (241, is still IR divergent and therefore a single quark 
is still absent from the spectrum. This conclusion sup¬ 
ports recent hadron spectroscopy studies, see Refs. Hi¬ 
lls, where the symmetries of hadrons are analyzed un¬ 
der low-mode truncation. The fact that confinement is 
still intact can be judged from the exponential decay sig¬ 
nal of hadron correlators. Here we have found a more 
physical, though gauge dependent, explanation for this 
phenomenon. 


VI. SUMMARY AND CONCLUSIONS 

A first detailed analysis of the overlap quark propag¬ 
ator in Coulomb gauge has been presented. The use 
of chirally symmetric lattice fermions makes the results 
highly predictive, since no improvement techniques have 
to be implemented, which could spoil the result. The 
drawback is the immense computational cost. To reach 
the small momentum region a large lattice has to be 
chosen, which is in our case a 20 4 lattice with a = 0.2 fm. 

We observe a clear indication of dynamical mass gener¬ 
ation in the quark propagator dressing functions. Scalar 
and vector dressing functions and therefore the dynam¬ 
ical quark mass start to increase in the momentum re¬ 
gion around 1 — 2 GeV. The constituent quark mass 
in the chiral limit is around 300 MeV, which is expec¬ 
ted to change only slightly for dynamical configurations. 
Moreover, we see hints that scalar and vector dressing 
functions, and therefore the energy dispersion, are di¬ 
verging in the infrared limit. However, such a state¬ 
ment can only be confirmed by taking the continuum 
limit. Moreover, for a suitable parameterization of the 
dynamical mass function, put forward in Ref. j4], the 
continuum limit is of great importance. Such a paramet¬ 
erization would be very helpful for continuum Coulomb 
gauge studies. 

We also studied the quark propagator after applying 
three levels of stout smearing to the gauge field configur¬ 
ations. Although the dressing functions show the same 
qualitative behavior as in the unsmeared case, they reach 
different small momentum values. The interesting obser¬ 
vation is that the dynamical mass function is not affected 
by the smearing procedure and is only shifted to slightly 
larger values. 
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In the second part of the paper we explored the in¬ 
terrelation between confinement and chiral symmetry 
breaking. Coulomb gauge is very appealing for this pur¬ 
pose, since an energy dispersion can be formulated, which 
relates the confinement properties of quarks to the di¬ 
vergence of the dressing functions. We created a scen¬ 
ario where the low modes of the Dirac operator are re¬ 
moved from the valence quark propagators and thereby 
the chiral condensate is suppressed. We found that, in 
contrast to the scalar dressing function B(p) : the vector 
dressing function A s (p) stays intact in the chiral limit 
and even seems to diverge stronger (again, for a conclus¬ 
ive result the continuum behavior has to be explored). 
Therefore the quark dispersion relation ui(p) can still 
be infrared divergent. It is the vector dressing function 
which is responsible that a single quark is still confined in 
the hadron. We stress that to our knowledge such a dir¬ 
ect connection between confinement and chiral symmetry 
breaking is only possible in Coulomb gauge. 

Moreover, the relation between the gauge dependent 
dynamical mass function and the gauge invariant chiral 
condensate has been given: the vanishing of the former 
implies the vanishing of the latter. 


The Coulomb gauge dispersion relation for confined 
quarks offers the possibility to study several interest¬ 
ing questions regarding the relation between confinement 
and chiral symmetry breaking, like, e.g., how do the 
dressing functions behave if instantons are removed from 
the quark propagator, or in the case of adjoint quarks? 
What happens at finite temperature? We expect that 
Coulomb gauge could shed light on these questions and 
help to understand the interrelation of confinement and 
chiral symmetry breaking. 
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